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ABSTRACT 

We present the results of an extensive observational study of the active star- forming complex W51 
that was observed in the J = 2 — 1 transition of the ^^CO and ^^CO molecules over a 1?25 x 1?00 
region with the University of Arizona Heinrich Hertz Submillimeter Telescope. We use a statistical 
equilibrium code to estimate physical properties of the molecular gas. We compare the molecular cloud 
morphology with the distribution of infrared (IR) and radio continuum sources, and find associations 
between molecular clouds and young stellar objects (YSOs) listed in Spitzer IR catalogs. The ratios 
of CO lines associated with H II regions are different from the ratios outside the active star- forming 
regions. We present evidence of star formation triggered by the expansion of the H II regions and by 
cloud-cloud collisions. We estimate that about 1% of the cloud mass is currently in YSOs. 
Subject headings: H II regions — stars: formation — ISM: individual (W51) — infrared: ISM — 
radio lines: ISM 
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1. INTRODUCTION 

Establishing the properties of molecular clouds is 
essential in understanding the formation and evolu- 
tion of star- forming regions. Newly formed mas- 
sive stars can affect the parental molecular clouds 
through ionization, heating, and expansion of the H II 
regions, stellar winds, and supernova- driven shocks 
(iBiedng et al."2QQ9': 'Beuther et al .1 120071: 1 Clark fc Ported 
20041: iGr itschneder et al. 2009; iPovich et all 120091 : 
Iborn et al.ll200 2). These mechanisms are either com- 
pressing or dispersing the surrounding clouds. For a bet- 
ter understanding of the feedback process in the interstel- 
lar medium (ISM), detailed observations of the molecular 
clouds and identification of embedded young stellar ob- 
jects (YSOs) are required. 

W51 is one of the most luminous star- forming re- 
gions in the first quadrant of the Galactic plane. 
The high lumi nosity c omes fr om a large number of 
0-type stars (jBiednd 119751 : iMehringed Il994l| that 
are within the m olecular cloud (JMufson fc Lisztl 119791 : 
iCarpenter fc Sand ers 199 8|). In th e 21 cm radio 
continuum map (Koo fc MoonI 119971 ). W51 can be 
divided into three regions: W51A (G49.4-0.3 and 
G49.5-0.4), W51B (G48.9-0.3, G49.0-0.3, G49.1-0.4, 
G49.10-0.27, G49.17-0.21,and G49.2-0 .4), and W51C 



(iKundu fc Velusamvlll96^lBiegind[T975l : lMufson fc Lisztl 
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While W51C shows a 



ri979l : iKoo fc Moor 
nonthermal con tinuum spectrum indicating a super- 
nova remnant ( Subrahmanvan fc GossI |1995), W51A 
and W51B have thermal spectra indicating compact 
H II regions. The ionized gas in W51A is dis- 
tributed in a broad velocity rang e of 40 to 80 km s~^ 
([Crampton et al. 1978; Mehringer 19941 : iPankonin et al] 
[19791 : I Wilson et al.l |1970). The average systemic ve- 
locity of the ionized gas in W51B is ~ 66 km s~^ 
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(iDownes et al.lll980l : [Pankonin et allll979l : IWilson et al.l 
Il970l ). In the Massachusetts-Stony Brook Galactic plane 
CO survey map cov ering the first Galactic quadrant 
(jSanders et al.l Il986f ) the W51 region is distinguished 
by bright CO emission extending over a 1° x 1° area 
centered on (/, b) ^ (49?5, — 0?2) at the velocity range 
from 4 5 to 75 km s~-^ (see ICar penter fc Sander si 1 19981 . 
Fig.l). ICarpenter fc SandersI ([1998) studied this region 
in detail in the ^^CO and ^^CO J = 1 - emission. 
They roughly divided the molecular clouds associated 
with the H II region complex into two groups: the 
W51 giant molecular cloud (CMC) and the 68 km s~^ 
cloud. They calculated a virial mass of 1.2 x 10^ Mq 
in the W51 CMC and 1.9 x 10^ M© in the 68 km 
s~^ cloud. Because the maxi mum velocity p ermitted 
by the Galactic rotation (Bra nd fc Blitzl Il993f ) toward 
the tangential point (/ = 49°) of the Sagittarius spiral 
arm is ~ 60 km s~^, the 68 km s~^ c loud is assoc i- 
ated with a high-velocity ( HV) stream ()Burtonlll970f ). 
ICarpenter fc SandersI (|1998l ) speculated that the massive 
star formation in W51 results from a collision between 
the W51 CMC and the 68 km s~^ cl oud whi ch is related 
to a s piral density wave. iKool (11999) and Kumar et al.l 
(|2004l ) also argued that enhanced star formation in W51 
might be caused by the spiral density wave. 

In this paper we present new observations of W51 in 
the J = 2 - 1 transition of ^^CO and ^^CO lines. Our 
on-the-fiy (OTF) maps are fully sampled in two lines si- 
multaneously with an angular resolution of 32^', and with 
rms noise per velocity channel of ~ 0.1 K in antenna tem- 
perature. Our sensitive, fully sampled mapping allows us 
to estimate mass more accurately and to delineate more 
detailed structure of the molecular gas. The purpose of 
this study is to find the correlation between the YSOs 
identified from the Galactic Legacy Infrared Mid-Plane 
Survey Extraordinaire (GLIMPSE I) Catalogs and the 
associated molecular clouds. We can approach the star 
formation history in this massive star-forming complex 
by analyzing CO data. In Section [2] we describe the ob- 
servations and data processing. In Section [3] we present 
our analysis of the spectral features and describe indi- 
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vidual regions based on the CO, radio continuum, and 
Spitzer data. In Section |4] we discuss the spatial distri- 
bution of the H II regions and estimate the star formation 
efficiency. In Section [5] we summarize our main conclu- 
sions. 

2. OBSERVATIONS AND PROCESSING 
2.1. 12(70 and ^^CO J = 2-1 and ^^CO J = 1-0 

The W51 GMC complex was mapped in the J = 2 — 1 
transition of ^^CO and ^^CO with the 10 meter Hein- 
rich Hertz Telescope (HHT) on Mount Graham, Ari- 
zona. The first observation was in 2005 November. A 
1?25 X 1?00 mapping centered on (/, b)= (49?375, -0?200) 
was completed in 2008 February. The whole map con- 
sists of a total of 15 sub-fields. Each of the 15' x 20' 
sub-fields was mapped with OTF scanning in right as- 
cension at 10'' per second with a row spacing of 10" in 
declination. Prior to 2006 December, we used a single- 
polarization double sideband SIS mixer receiver with sys- 
tem temperatures of typically 350 - 500 K, depending on 
weather and source elevation. Since 2006 December, the 
observations were made using the 1.3 mm ALMA dual- 
polarization sideband-separating receiver with a 4-6 GHz 
intermediate frequency band. The receiver was tuned to 
the 12 CO J = 2 - 1 line at 230.538 GHz in the upper 
sideband and the ^^CO J = 2 - 1 line at 220.399 GHz in 
the lower sideband. Typical system temperatures were in 
the range 180 - 300 K. The spectrometers, one for each 
of the two polarizations and the two sidebands, were fil- 
ter banks with 256 channels of 1 MHz width and sepa- 
ration (corresponding to a velocity resolution of 1.3 km 
s~^ for the ^^CO J = 2 — 1 line). For some subfields of 
the ^^CO J = 2 — 1 mapping, the spectra were measured 
with an acousto-optic spectrometer (AOS) having a fre- 
quency resolution of 1 MHz, and a total bandwidth of 1 
GHz. The main beam efficiency, measured from planets, 
was 0.85 with an estimated uncertainty of (+0.05,-0.10). 

The raw data were processed with the CLASS reduc- 
tion package (from the University of Grenoble Astro- 
physics Group). The intensity scales for the two po- 
larizations were determined from observations of W51D 
made just before the OTF maps. System temperatures 
were chec ked by the standard am bient temperature load 
method (IKutner fc Ulichl Il981f ) after every other row 
of the map grid. Further analysis was don e with the 
MIRIAD software package (|Sault et al.lll995[ ). Although 
the data observed early in the program with a differ- 
ent receiver were noisier than the data obtained using 
the ALMA receiver, data quality was improved after 
combining all of the data using a {1/Tsys)'^ weighting. 
At the final step, the data were convolved with a 16" 
(FWHM) Gaussian beam, giving a 36" effective resolu- 
tion from the original angular resolution of the telescope, 
32" (FWHM). The two polarizations were averaged for 
each sideband, yielding images with rms noise per ve- 
locity channel of 0.16 and 0.07 K in T^ for the ^^CO 
J = 2 — 1 and ^^C0J = 2 — 1 respectively. The average 
rms noise level of the ^^CO map is larger than ^^CO be- 
cause several sub-fields of ^^CO map were observed just 
once with the noisier receiver. The full set of HHT data 
will be made available in a separate paper (Bieging et al. 
2010, in preparation). 

We obtained ^^C0J = 1 — data for the same region 



from the Galactic Ring Survey (GRS) ^ (jJackson et al.l 
l2006f ). Their spectral resolution is 0.2 km s~^, and the 
angular resolution and sampling are 46" and 22", re- 
spectively. We transform antenna temperature of ^^CO 
J = 1 — of GRS to main beam temperature, dividing 
by a main beam efficiency of 0.48. 

2.2. Spitzer Data 

The Galactic Legacy I nfrared Mid- Plane Su rvey Ex- 
traordinaire (GLIMPSE L iBeniaminet al.ll2003[ ) covered 
the Galactic plane (10° < |/| < 65°, \b\ < 1°) with the 
four mid-IR ban ds of the Infrared Array Camera (IRAC; 
iFazio et al.l 1200 4) on the Spitzer Space Telescope. Each 
IRAC band contains different spectral features : Band 1 
(3.6 /im) shows mainly continuum emission from stars. 
Band 2 (4.5 /im) traces H2 rotational transitions arising 
in the shocked gas associated with outflows, and emis- 
sion in Bands 3 (5.8 /im) and 4 (8.0 /im) is dominated 
by polycyclic aromatic hydrocarbon (PAH) features. For 
this study, we have retrieved images of the 1?25 x 1?00 
region of W51 centered on (/,6) = (49?375, -0?200) by 
combining two mosaic images which have a high resolu- 
tion of 1'.'2. For the same region, we used the GLIMPSE 
I Catalog to extract point sources that are detected at 
least twice in one band with an S/N > 5. The GLIMPSE 
I Catalog also tabulates JHKg flux densities from the 
Tw o Micron All Sky Surv ey (2MASS) point source cata- 
log (jSkrutskie et al.ll2006f ) for ah GLIMPSE sources with 
2MASS identi fications. 

MIPSGAL (Carey et al."2005') is a legacy program cov- 
ering the inner Galactic plane, 10° < |/| < 65° for 
\h\ < 1°, at 24 and 70 /im wit h the Multiband I maging 
Photometer for Spitzer fMIPS: iRieke et al.l[20Ql ) Space 
Telescope. The resolution of the 24 /im mosaics from 
the MIPSGAL survey is 2'.'4. We extracted 24 /im point- 
sources with F/SF > 7, then bandmerged the 24 /im 
sources with the GLIMPSE Catalog sources using a 2'.'0 
correlation radius. The final source list contains data 
from ah 8 bands: 2MASS {JHKs), GLIMPSE (IRAC 1, 
2, 3, 4 bands), and MIPSGAL 24 /im. A total of 104,582 
sources within the target region were select ed from the 
GLIMPSE I Catalog. In an earlier paper, 'Kang e t al.l 
(2009b) examined this set of point sources to identify 
and classify YSOs near the W51 H II region complex us- 
ing spectral energy distribution (SED) fits. All 8 bands 
are used for SED fitting as long as detected. 

3. RESULTS 
3.1. Velocity Structure 

Toward the W51 H II region complex, most ^^CO emis- 
sion comes from velocity i ntervals of Vlsr == 0-2 5 km 
s-^ and 35-75 km s"^ (Car penter fc Sander si [T998h . We 
do not consider the former velocity component, which 
originates from nearby molecular clouds. We focus on 
the later velocity component, which is extended from 30 
to 85 km s~^ in our sensitive and small grid mapping, 
associated with the W51 H II region complex. 

Figure [H represents the ^^CO J = 2 — 1 intensity map 
integrated over the velocity range between 30 and 85 km 
s~^. Most of the bright ^^CO J = 2 — 1 emission peaks 
are coincident with bright radio continuum sources. The 

^ http://www.bu.edu/galacticring/ 
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Fig. 1. — Distribution of the ^-^CO J = 2 — 1 emission in the 
W51 H II region complex. The ^"^CO J = 2 — 1 emission is inte- 
grated from 30 to 85 km s~^. The contour levels are 1, 5, 10, 30, 
50, 70, and 90% of the peak value (542 K km s~^) in the map. 
Star symbols mark the thermal radio continuum sources listed by 
iGoss & Shaver (1970). The names of the H II regions are G49.5- 
0.4, G49.4-0.3, G49.2-0.3, G49.1-0.4, and G48.9-0.3, from east to 
west. 

W51 H II region complex shows a complicated structure 
in wide velocity ranges. Carpenter fc Sand ers (1998) di- 
vided molecular clouds associated with W51 into seven 
components at 44, 49, 53, 58, 60, 63, and 68 km s~^. 
In W51A, the 53 km s~^ component is associated with 
G49.4-0.3 and the 58, 60, 63 km s~^ components are as- 
sociated G49.5-0.4. The 68 km s~^ cloud is related to 
W51B. Besides the bright emission associated with radio 
continuum sources, there is ^^CO J = 2 — 1 emission ex- 
tended over the whole of the 1?25 x 1?00 region. More 
detailed structures of ^^CO and ^^CO J = 2 — 1 emission 
are shown in channel maps (Figures [2] and [3]). 

In the 43 and 46 km s~^ channel maps of Figure [21 
there are molecular clouds not associated with the ra- 
dio continuum sources. Some clouds are parallel to the 
Galactic plane at 6 = 0°, and the others in the western^ 
part of the map are above the Galactic plane. CO emis- 
sion in the W51A region appears in the first channel, 
43 km s~^. Especially at the position of the brightest 
radio continuum source, G49.5-0.4, the CO emission is 
distributed through the whole velocity range from 46 to 
76 km s~^. In contrast, the continuum sources of W51B 
are coincident with the filamentary structure with veloc- 
ity > 6 4 km s~^, which we d iscuss below. As pointed 
out by [Carpenter fc SandersI (|1998f ), the W51 molecu- 
lar clouds with velocities of 58, 63 km s~^ are extended 
through a large area of the map, and the southern edge 
of the diffuse emission seen in the 64 km s~^ channel map 
is seen against the filamentary structure of 67 km s~^. 
Two holes appear in the northern (/,6 = 49?50, — 0?30) 
and eastern (/, b = 49?60, -0?35) parts of G49. 5-0.4 from 
49 to 64 km s^^ (Figure [2j). 

Figure m shows the total CO distribution around W51 
in three velocity ranges, 30 — 55, 56 — 65, and 66 — 85 

^ Note: In this paper, directions are with reference to galactic 
coordinates, i.e., "north" means toward increasing latitude and 
"east" means toward increasing longitude. 
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Fig. 2. — Channel maps of the ■'^^CO J = 2 — 1 emission averaged 
in 3 km s"-*^ intervals over the velocity range from 43 to 76 km 
s"-*^. The values printed in the top left corner of each panel denote 
the center of the velocity intervals. The gray scale in each panel 
is a logarithmic stretch from log 1.0 K (white) to log 10.0 K in 
Tmb (black). Open circle symbols mark the thermal radio sources, 
and cross symbol represents the non-thermal radio source listed by 
Goss &: Shaver (J970). 
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Fig. 3. Same as Fig. ^ for ^^CO J = 2-1. The gray scale 
in each panel is a logarithmic stretch from log 0.3 K (white) to log 
4.0 K in Tmb (black). 

km s~^. The molecular clouds associated with each ra- 
dio continuum source show a range of different velocity 
structure. The white boxes outline the regions that we 
have analyzed to estimate the physical quantities from 
the i^CO and ^^CO J = 2 - 1 data. 
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TABLE 1 

Physical parameters^ 



Name^ 




Selected 


area 




Velocity range^ 


Mass^ 


Ay- 


^^i^2-l/l-0^ 


13/12^^2-1^ 




^minV ) 


^max(°) b 


mini, ) 


6max(°) 


(km s-i) 


(lO^Mo) 


(mag) 






G49.5-0.4 


49.42 


49.60 


-0.45 


-0.29 


49 - 72 


8.7 


72 


0.87 


0.31 


G49.4-0.3 


49.32 


49.44 


-0.36 


-0.23 


47-55 


3.7 


26 


1.01 


0.26 


G49.2-0.3 


49.13 


49.29 


-0.42 


-0.27 


62-73 


3.8 


21 


0.78 


0.25 


G49.1-0.4 


49.04 


49.13 


-0.42 


-0.31 


56-76 


1.6 


20 


0.98 


0.23 


G49.17-0.21 


49.13 


49.21 


-0.24 


-0.18 


58-69 


0.9 


20 


0.78 


0.37 


G49.10-0.27 


49.06 


49.14 


-0.31 


-0.24 


58-70 


1.3 


22 


0.83 


0.25 


G49.0-0.3 


48.96 


49.03 


-0.34 


-0.26 


62-75 


1.5 


26 


1.14 


0.21 


G48.9-0.3 


48.85 


48.96 


-0.32 


-0.24 


63-74 


1.9 


26 


1.16 


0.20 


Shell 


49.60 


49.77 


-0.24 


-0.06 


55 - 68 


2.9 


15 


0.43 


0.12 


Cluster 


49.30 


49.55 


-0.13 


0.12 


59-69 


6.0 


15 


0.41 


0.13 


IRDCl 


48.75 


48.88 


0.06 


0.17 


47-56 


1.7 


14 


0.54 


0.17 


IRDC2 


48.80 


48.93 


0.18 


0.28 


42-49 


1.1 


9 


0.52 


0.17 



We assumed a distance of 6 kpc from the Sun to derive the physical parameters. 

See Figure U 

Velocity range considered in the estimation of mass and line ratios 

Mass derived using an LVG escape probability radiative transfer model 

Peak Av derived from the CO column density 

Mean ^^COJ = 2-1/J = 1-0 intensity ratio 

Mean ^^CO / ^^CO J = 2 - 1 intensity ratio 
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Fig. 4. — Color image of the W51 region composed of -"^^CO 
J = 2 — 1 intensity integrated over the velocity range of 30 — 55 km 
s"-*^ (blue), 56 — 65 km s"-*^ (green), and 66 — 85 km s"-*^ (red). Cross 
symb ols represent the compact radio continuum sources listed by 
IKool ([1997). White boxes are the area for estimating the physical 
parameters in Table [T] 

3.2. Line Ratios 

Figure [5] shows the average spectra and ratios for ah 
the regions marked with white boxes in Figure |H CO 
hue ratios are useful indicators of trends in hue opacity 
and gas column density (isotopologue ratios, e.g., ^^CO / 
^^CO J = 2 - 1 (= ^^/^2i?2-i)); or of molecular excita- 
tion (e.g., i^CO J = 2-l/J = l-0(= ^^i^2-i/i-o)). 
As seen in Figure O G49.5-0.4 has a wide velocity range, 
from 45 to 75 km s~^ showing a peak at 60 km s~^. The 
average ^^i?2-i/i-o value is 0.87 within the two vertical 
dotted lines. The selected velocity ranges, covering the 
strongest ^^CO emission, are associated with the H II 
regions or YSOs. The ^^/^^i?2-i value increases with ve- 
locity up to around 65 km s~^ and then decreases. The 
spectrum of G49.4-0.3 shows a peak at a lower velocity, 
51 km s~^. The spectrum in W51B has strong emis- 



sion between 55 and 75 km s ^, which is different from 
the molecular cloud associated with the W51A region. 
The ^^i^2-i/i-o values of G49.2-0.3, G49.1-0.4, G49.10- 
0.27, G49.0-0.3, and G48.9-0.3 in W51B systematically 
increase in the higher velocity ranges (Figure [5]). The 
ratio plots in Figure [5] generally show an anti-correlation 
between ^^^^^R2-i and ^^i?2-i/i-o. In Table[I]we list the 

average ^^i?2-i/i-o and ^^/^^i?2-i values derived from 
each region. 



The average R2 



-1/1-0 



ratios associated with all H II 



regions are greater than 0.7, which implies that these 
clouds are relatively warm and dense. In contrast, 
the ratios of ^^i?2- 1/1-0 foi" the molecular clouds out- 
side the active star-forming region, where there are no 
known continuum sources, are smaller than 0.7 (see Shell, 
Cluster, IRDCl, and IRDC2 in Figure E]). Molecu- 
lar clouds with low ratios were often observed in dark 
clouds and outer env elopes of giant molecular clouds 
(jSakamoto et al.lll997[ ). 

3.3. Physical Properties 

We use a non-LTE statistical equilibrium treatment of 
the CO molecular excitation to study physical proper- 
ties of the molecular cloud. We apply the escape prob- 
ability radiat i ve tra nsfer and photodissociation model of 
iKulesa et al] (|2QQ5l ) to our ^^CO and ^^CO J = 2 - 1 
data. We calculate grids of CO level populations for 
wide ranges of volume densities (10^ — 10^ cm~^) and 
temperatures (5 — 300 K) assuming detailed balance and 
steady state. From these model grids, the total CO col- 
umn density at each observed pixel is computed from 
the peak temperature, line widths and integrated inten- 
sities of the observed CO lines. The velocity range for 
integrating CO lines was determined from the range cov- 
ering the strongest ^^CO J = 2 — 1 emission on the line 
of sight (Figure [5j). Assuming that the CO heating is 
dominated by photon processes (e.g., the photoelectric 
heating of dust), a coarse estimate of the incident radi- 
ation field can be made for each point in the map. The 
photodissociation model is then applied to estimate a to- 
tal hydrogen column density from the CO data. This 
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Fig. 5. — Average spectra of ^^CO J = 2 - 1 
■'^^i^2-i/i-o (r^d crosses) and ^^ I ^'^ R^-x (green 
vertical dotted lines show the velocity range for 
is smaller than marker. 



(thick hne), ^^CO J = 2 - 1 (thin hne), and ^^CO J = 1 - (dashed hne) and ratios of 
dots) for the boxed areas in Figure |4] The horizontal dotted line is at 0.7 (see text). Two 
deriving mass and the average line ratio of clouds, listed in Tabled A typical uncertainty 
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Fig. 6. — Top: Color image of the W51A region composed of 
IRAC 5.8 /im (red), 4.5 /im (green), and 3.6 /im (blue). White 
boxes show the area for analyzing the labeled radio continuum 
source in detail. Middle: Contour map of -"^^CO J = 2 — 1 intensity 
integrated from 45 to 65 km s"-*^. The contour levels are 40, 80, 
120, 160, 200, 240, 280, 320, 360, and 400 K km s-^. The oblique 
line shows the cut for the PV maps presented in Figure [7] Bottom: 
Squares are compact radio continuum sources listed by Mehringer 
([l994i ). YSO candidates (open circles) are marked in red for Stage 
O/I, yellow for Stage II, and green for ambiguous sources. 
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Fig. 7. ^^qq J = 2 - 1 (top), ^^co J = 2 - 1 (middle) and 
12/I3j^2_i (bottom) position- velocity map through the W51 A 
region from southeast (A) to northwest (A^) marked at the second 
panel of Figure [6] The temperature scale is indicated on the 
right (K in Tmb)- Contour levels of ^^ CO J = 2 - 1 and ^^CO 
J = 2 — 1 maps are 10, 20, 40, and 80% of peak value in the map. 
Contour levels of ^^^^^R2-i map are 1.5, 3.0, 5.0, and 7.0. Star 
symbols represent the position of the radio continuum sources (d 
and b of G49.5-0.4 and b, e, and d of G49.4-0.3). The velocities of 
the radio continuum sourc es come from the radio recom bination 
line observations (|van Gorkom et al.lll980l : IPankonin et al.iiil979ii V 



calculation is based on the CO and H2 photo di ssoci- 
ation treatments of | van Dishoeck fc Blac3 (|1988f ) and 
iBlack fc van Dishoeckl ()1987[ ). respectively, using a to- 



tal interstellar carbo n abundance of C/H = 2.4 x 10~^ 
([Cardelli et al.ilT996[ ). Because the CO abundance is a 
strong function of column density, blind application of a 
uniform CO "dark cloud" abundance (10""^) can lead to 
a gross underestimate of total hydrogen column density 
and gas mass. The visual extinction, Ay, is calculated 
from N(H) usingj; N (H)/Av = 1.9 x 10^^ cm-^ mag"^ 
(|Bohlinet al.lll978f ). The distance to W5 1 derived in pre- 
vious studies vary from 2.0 to 8.3 kpc. iFigueredo et al.l 
(|2QQ8f ) derived 2.0±0.3 kpc by spectroscopic parallaxes 
of 0-type stars. The kinematic distance is 5.5 kpc using 
radio recombination lines (.Russei l 2003.) . The values de- 
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Fig. 8. — Top: Color image of the W51B region composed of 
IRAC 5.8 /im (red), 4.5 /im (green), and 3.6 /im (blue). White 
boxes show the area for analyzing the labeled radio continuum 
source in detail. Middle: Contour map of -"^^CO J = 2 — 1 intensity 
integrated from 55 to 75 km s"-*^. The contour levels are 30, 90, 
150, 210, and 270 K km s"-*^. The oblique line shows the cut for 
the posit ion- velocity maps presented in Figure [9] Bottom: Squares 
are compact radio continuum sources listed by Koo (1997). YSO 
candidates (open circles) are marked in red for Stage O/I, yellow 
for Stage II, and green for ambiguous sources. 

rived from water maser proper motion measurements are 
6.1±1.3 kpc (Imai et al. 2002), 7±1.5 kpc (Genzel et alJ 
|^81D, an d 8.3±2.5 kpc (Schneps et al.lll98 1). Recently, 
IXu et al] (|2009 ) found the distance of 5.ll^;4 kpc using 
trigonometric parallax measurements. In this paper, we 
adopt a distance of 6 kpc and use that value to calculate 
the gas masses of the selected regions shown in Figure |H 
In Table [H we list the size of the regions, velocity range, 
masses, extinctions, and the mean ratios of the lines. 

3.4. Large Scale Structure 

3.4.1. W51A 

W51A includes the radio continuum sources G49.5- 
0.4 and G49.4-0.3 (Figure [6]). 049. 5-0.4 is the bright- 
est component in radio continuum in the W51 re- 
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Fig. 9. ^^qq J = 2 - 1 (top), i^co J = 2 - 1 {middle) and 
12/I3j^2_i (bottom) position- velocity map through the radio con- 
tinuum sources G49.2-0.3, G49.0-0.3, and G48. 9-0.3 from B to B^ 
in Figure [8] The temperature scale is indicated on the right (K in 
Tmb)- Contour levels of the ^^ CO J = 2 - 1 and ^^ CO J = 2 - 1 
maps are 10, 20, 40, and 80% of the peak value in the map. Contour 
levels of the ^2/13^2-1 map are 1.5, 3.0, 5.0, and 7.0. Star sym- 
bols represent the radio continuum sources. The RRL velocities of 
radio continuum sourc es are the mean velocities of the ionized gas 
(IPankonin et al.llT979h. 



gion. It contains the infrared so urces W51 IRSl and 
IRS2 (Wv nn- Williams et al.l[T97 4^. and the H2O masers 
W51 MAIN W51 N, and W51SlGenzel & Downes 19^^. 
iMehringerl (jl994[ ) found that G49. 5-0.4 consists of 16 dis- 
crete components using the VLA radio continuum obser- 
vations at 3.6, 6, and 20 cm. G49.4-0.3 located to the 
west of G49.5-0.4 i s associated with 6 radio continuum 
sources (|Mehringer 1994). 

Figure [6] shows the distribution of mid-infrared and 
CO emissio n in W51A . The radio continuum sources 
detected by IMehringerl (|1994f ) are marked with squares 
and YSO candidates are marked with colored open cir- 
cle s. YSO c andidates were identified and classified 
by Kang et al. (2009b) using the SED fitting tools of 
Robitaille et aL (2007) . Each color represents the evo- 
lutionary stage of the YSO candidates: red for Stage 
I/O, yellow for Stage II, and green for ambiguous sources. 
Most YSO candidates in the W51A region are Stage 
I/O (red) or ambiguous sources (green). The bright 
IRAC 4.5 jam emission (green) is seen at the position 
of the radio continuum sources. This 4.5 jam emission 
might be explained as l ine emis sion fro m H2 in shocked 
molecular ^as fl Cvgano wski et al . 2008; Davis et al.ll2007l : 
IShepherd et~al.ii2007i : iSmith et al.i 12006). 

Figure [7] shows ^^CO J = 2 - 1, ^^CO J = 2 - 1, and 
the ratio of ^^CO J = 2-l to ^^CO J = 2-l (= 
12/13^^_^^ in the form of position- velocity (PV) maps of 
G49.5-0.4 and G49.4-0.3 along the line shown in Figure 
[6] (middle). The PV maps illustrate the kinematic rela- 
tionship of each radio continuum source and associated 
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molecular cloud in W51A. The radio continuum sources 
of G49.5-0.4 appear to be associated with the molecular 
cloud at 58 km s~^, and the G49.4-0.3 sources with a 53 
km s~^ molecular component. Figure [71 also shows the 
dumpiness of the HV stream as well as the spatial re- 
lationship between the radio continuum sources and the 
molecular clouds. The HV stream with velocity 68 km 
s~^ is strong in front of the G49.5-0.4 region, while the 
HV stream near G49.4-0.3 is weak. 



3.4.2. W51B 

W51B consists of six radio continuum sources: G49.2- 
0.3, G49.17-0.21, G49.10-0.27, G49.1-0.4, G49.0-0.3 and 
G48.9-0.3. In Fi gure [ H I we mark the radio continuum 
sources listed bv lKool() 19971 ) on the three-color IRAC 
composite image. The molecular cloud toward the W51B 
region is distributed within the velocity range from 55 to 
75 km s-^ Figure [8] (middle) shows the ^^CO J = 2 - 1 
intensity integrated between 55 and 75 km s~^, which is 
coincident with the bright IRAC emission. Most strong 
molecular clouds appear around the radio continuum 
sources. Eighty percent of the YSO candidates (lower 
panel) in W51B are early stage objects (Stage O/I or 
ambiguous). 

The H II regions of W51B are distributed parallel to 
the Galactic plane. Molecular clouds are distributed 
along 6 ~ — 0?3. To analyze the detailed structure of 
each H II region and associated molecular cloud, we di- 
vide the W51B region into six groups: G49.2-0.3, G49.17- 
0.21, G49.10-0.27, G49.1-0.4, G49.0-0.3 and G48.9-0.3. 

Figure [9] shows the PV maps of ^^cq J = 2 - 1, ^^CO 
J = 2 — 1 and ^^/^^i?2-i, through the three H II regions 
in the W51B region from B to B' in Figure [8] (mid- 
dle). Most CO emission appears between 55 and 75 km 
s~-^. The molecular clouds associated with G48.9-0.3 and 
G49.0-0.3 are strong near 68 km s-^ The ^2/^^i?2-i ra- 
tio around G49.2-0.3 is smaller than 1.5 (red contours) 
at the offset position of — 7.'5, which implies that the 
optical depth of the molecular gas is large at that posi- 
tion. Several YSOs are identified to the east of G49.2- 
0.3 {l,h = 49?26,-0?34) on Figure [S In Figure [9] (bot- 
tom), the molecular cloud associated with these YSOs, 
also shows a low ratio (^^/^^i?2-i < 1.5), like the molec- 
ular cloud associated with G49.2-0.3. 

3.5. Individual Regions 

3.5.1. G49.5-O.4 

Figure [10] shows the complicated structure of G49.5- 
0.4 in various wavelengths. Although the IRAC emission 
and the ^^CO J = 2 — 1 integrated intensity show sim- 
ilar structure, the peak position of the ^^CO J = 2 — 1 
intensity coincides with a dark part of the IRAC com- 
posite image near the southeast region of component b. 
Some radio sources are located in the brightest part of the 
molecular cloud. Two bright radio continuum sources, 
W51d and W51e, coincide with the infrared sources IRS 
2 and IRS 1 and are adjacent to two peaks of the in- 
tegrated intensity map of the ^^CO J = 2 — 1 emission. 
Other sources (a, b, f, g, h, and i) are located near the 
edge of the molecular cloud. These radio sources are co- 
incident with the bright region in the IRAC 4.5 jam band 
(green) probably indicative of shocked H2. 
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Fig. 10. — The ^^CO J = 2 - 1 intensity integrated over the ve- 
locity range from 45 to 75 km s"-*^ overlaid on the color image of the 
G49. 5-0.4 region composed of IRAC 5.8 //m (red), 4.5 /im (green), 
and 3.6 /im (blue). Contour levels are 20, 40, 60, 80, 100, 120, 140, 
160, 180, 200, 220, and 240 K km s-^ Squares are the ZAMS OB 
stars listed by Mehringer (1994). YSO candidates (open circles) 
are marked in red for Stage O/I, yellow for Stage II, and green for 
ambiguous sources. The dash-dotted rectangle shows the field in 
Figure fTTl 
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Fig. 11. Tile map of the ^^CO J = 2 - 1 (thick solid line), 
13 CO J = 2 - 1 (thin solid fine), and ^^ CO J = 1 - (dash-dotted 
line) spectra of G49. 5-0.4 component b seen in the dash-dotted box 
in Figure [To] Spacing of the mapping grid is 22^^ 

The spatial distribution of the radio continuum sources 
and associated molecular clouds can be constructed by 
comparing the velocity of ionized gas, CO emission, and 
absorption lines. The average velocity of the ionized gas 
is 59 km s~^, whic h was measured usin g the H109a re- 
combination line by I Wilson et al.l (jTOTOl). CO emission is 
brightest at Vlsr = 61.0 km s~^ (Figure[5]). The ionized 
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Fig. 12. — Spatial and posit ion- velocity diagrams of component 
b of G49. 5-0.4 in ^^CO J = 2 - 1. (a) is channel map at 64 km 
s"-*^, which is RRL velocity. The star symbol represents component 
b of G49. 5-0.4. (b) and (c) are position- velocity maps in latitude 
and longitude across the solid lines in (a). The ■'^^CO contour levels 
are 1.6, 3.2, 6.4, 12.8, and 25.6 K in Tmb- (d) is composite image 
of IRAC 5.8 /im (red), 4.5 /im (green), and 3.6 ^ m (blue). The 
contours show the 21 cm radio continuum emission (|Koo &: MoonI 
fT997i) . The contour levels are 0.015, 0.05, 0.1, 0.15, 0.3, 0.6, 0.9, 
1.2, 1.5, and 1.8 Jy beam"!. 
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Fig. 13. Same as in Figure fT2l for I^CO J = 2 - 1. The ^^CO 
contour levels are 0.5, 1.0, 2.0, 4.0, 8.0, and 16.0 K in Tmb- Con- 
tours on (d) show the -"^^CO J = 2 — 1 intensity integrated from 42 
to 75 km s-i. The contour levels are 15, 45, 75, 105, 135, 165, 195, 
225, 255, and 285 K km s-^. 
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Fig. 14. — Same as in Figure[T2] for ratio of ^^/^^i?2-i- Contour 
levels are 1.5, 3.0, 5.0, and 7.0. YSO candidates (open circles) on 
(d) are marked in red for Stage O/I, yellow for Stage II, and green 
for ambiguous sources. 



and molecular gas have similar velocities, which implies 
that the molecular cloud at this velocity is associated 
with the radio continuum source. Toward G49.5-0.4, the 
central velocities of the ionized gas vary from 48 to 77 
kms~^ (M ehringerl ' 19941 . H2CO a bsorption has been de - 
tected between 52 and 73 km s~^ (jArnal fc Gossl[l985f ). 
H I 21 cm absorption has been observed at 50, 62, and 



69 km s 



a psorp tK 

(HEBE 



1997i) . H2CO absorption between 66 



and 70 km s~^ and H I absorption at 69 km s~^ implies 
that the HV stream is in front of G49.5-0.4. Regard- 
ing the spatial location of each radio continu um source 
in the G49. 5-0.4 region, lArnal fc GossI (J1985D suggested 
that components a, b, and c are closer to the observer 
than components d and e by c ompa r ing th e Ha 109 RRL 
and H2CO absorption spectra. iKool (11997) reported that 
the H I cloud with 62 km s~^ velocity associated with 
G49.5-0.4, and components a, b, and e are either embed- 
ded in or behind this cloud while components f and g are 
located either in front of or the outside the cloud. Also, 
all radio continuum sources in the G49.5-0.4 region are 
spatially close because all H I spectra detected toward 
these sources have similar absorption profiles. Most of 
the bright radio components of G49.5-0.4 in the W51 
H II region complex must be directly associated with the 
molecular clouds at 60 km s~^. 

CO spectra in the G49.5-0.4 region show complex 
structure. The average ^^CO J = 2 — 1 spectrum of 
G49.5-0.4 in Figure [5] is divided into two components 
with velocities of 59 and 68 km s~^ at the emission 
peaks. The brighter component between 45 to 65 km 
s~^ is associated with the radio continuum sources in 
G49.5-0.4. Figure [TTI shows the CO spectra around com- 
ponent b of G49.5-0.4. There are three peaks in the ^^CO 
J = 2 — 1 emission at 52, 58, and 70 km s~^. The inten- 
sities of ^^CO J = 2 - 1 and J = 1 - at 66 km s-^ are 
greater than ^^CO J = 2 — 1 intensities, which implies 
self- absorption in the ^^CO J = 2 — 1 emission. These 
CO self-absorption features clearly appear as a low ratio 
of ^^/^^i?2-i on PV maps inside the first contour (< 1.5) 
o n Figure [71 (bot t om) . 

lArnal fc GossI (|1985l ) proposed that component b of 
G49.5-0.4 is an H II regions created by O stars formed 
in the collision between the HV stream and the molec- 
ular cloud at 61-63 km s~^. The RRL velocity of com- 
S )onent b has been detec ted at 63.8 km s~^ in H109a 
van Gorkom et al] |1980[ ) . Figures [T2HT41 show the spa- 
tial and kinematic structure around component b of 
G49.5-0.4 in ^^CO, ^^CO J = 2 - 1, and i^/is^^^-i. As 
seen in Figures [T2HHI component b lies between the 
main molecular cloud and the HV stream on the PV 
map. A clear difference between the PV maps of each 
CO isotopologue appears at 66 km s~^. The ^^CO 
J = 2 — 1 intensity is very weak at the velocity between 
the main molecular cloud and the HV stream. How- 
ever, the PV map of^^CO J = 2 — 1 does not show an 
apparent boundary at this velocity so that ^^/^^i?2-i be- 
comes small near 66 km s~^. A d eep abso rption fea- 
ture of H2CO spectra (|Fomalont fc W ehachewl ll973l ) and 
lower opacity molecular species such a s CS J = 3 — 2 
(Penzias et al. 1971) and HCN ( Snvd eTlFBuhl [T97ll ) 
have been detected at 66 km s~^, which imply the pres- 
ence of cold dense material in this region. In Figure 
[Hlfd) there are dark filamentary structures on the IRAC 
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Fig. 15. — Spatial and posit ion- velocity diagrams of component 
g of G49.5-0.4 in ^^CO J = 2 — 1. (a) is the channel map at 55 km 
s"-*^. The star symbol represents the component g. Contour levels 
are 1.5, 3, 6, and 12 K in Tmb- (b) and (c) are position- velocity 
maps in latitude and longitude across the solid lines in (a), (d) is 
the composite image of IRAC 5.8 /im (red), 4.5 /im (green), and 
3.6 /im ( blue). The contour s show the 21 cm radio continuum 
emission (|Koo fc Moonlll997l V The contour levels are 0.01, 0.03, 
0.05, 0.10, 0.15, and 0.18 Jy beam"!. Scale bar in the bottom left 
corner represents 1 pc at the distance of 6 kpc. 
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Fig. 16. — The ^^CO J = 2 - 1 intensity integrated over the ve- 
locity range from 45 to 55 km s"-*^ overlaid on the color image 
of the G49.4-0.3 region composed of IRAC 5.8 fim (red), 4.5 /im 
(green), and 3.6 /im (blue). Contour levels are 25, 50, 75, 100, 125, 
150, and 175 K km s"-*^. Squares are the ZAMS OB stars listed 
by Mehringer (1994). YSO candidates (open circles) are marked 
in red for Stage O/I, yellow for Stage II, and green for ambiguous 
sources. The oblique lin e sh ows the cut for the posit ion- velocity 
map presented in Figure fTTl 



composite image. Several YSOs in very early stages are 
detected in the dark filaments which coincide with the in- 
ner area with ^^/^^i?2-i of 1.5 (the first contour in Figure 

ma)). 

Component g of G49.5-0.4 was observed at 54 km 
s~^ in [S III] emission line (Goudis & Hippelein 1982). 
The average Ha velo city is 54.6 ± 13.7 km s~^ 
(|Crampton et al.l Il978f ). To investigate the velocity 
structure around component g, we present in Figure [15] 
the channel map at 55 km s~^, which is the ionized gas 
velocity, and the position velocity map across compo- 
nent g. A prominent molecular gas cavity exists around 
component g in Figure [TSTa). This cavity appears at 
a velocity of 55 km s~-^ of both (b^v) (b) and {l^v) (c) 
maps. Figure [TSTc) shows the velocity gradient from the 
expansion of molecular clouds around component g. In 
addition to the velocity structure, the green color in Fig- 
ure [TST d) caused by the shocked H2 emission is a strong 
evidence for the interaction of the ionizing star(s) and 
surrounding molecular cloud. 

3.5.2. G494-O.3 

G49.4-0.3 is located to the northwes t of G49. 5- 0.4 and 
is composed of six subcomponents (|Mehringerl [1994). 
Figure [16] shows the detailed structure of G49.4-0.3 in 
12 CO J = 2 - 1 and IRAC bands. The bright emission 
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Fig. 17. ^^qq J = 2-1 (left), ^^CO J = 2-1 (middle) and 
12/I3j^2_]^ (right) position- velocity map through the radio con- 
tinuum sources a, e, and f of G49.4-0.3 from southeast (C) to 
northwest (C^) of the map in Figure [16] The temperature scale 
is indicated on the top (K in Tmb)- Contour levels of the -"^^CO 
J = 2 - 1 and I^CO J = 2-1 maps are 10, 20, 40, and 80% of 
peak value in the map. Contour levels of the ^^/ ^^i?2-i map are 
1.5, 3.0, 5.0, and 7.0. Star symbols represent the position of the ra- 
dio continuum sour ces. The velocity of a and e is the mean velocity 
of the ionized gas (jPankonin et al.lll979l ). The velocity of compo- 
nent f comes from the velocity of the molecular cloud around the 
radio continuum source. 
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Fig. 18. — The ^^CO J = 2 - 1 intensity integrated over the 
velocity range from 57 to 76 km s"-*^ overlaid on the color image 
of the G49.2-0.3 region composed of IRAC 5.8 /xm (red), 4.5 /im 
(green), and 3.6 /im (blue). Thick contour levels are 7, 14, 21, 
28, 35, and 42 K km s"-*^. Thin white contours show 21 cm radio 
continuum map. Contour levels are 0.015, 0.03, 0.06, 0.1, 0.3, 
0.6, 0.9, 1.2, 1.5, and 1.8 Jy beam"! (Koo & Moon 1997). YSO 
candidates (open circles) are marked in red for Stage O/I, yellow 
for Stage II, and green for ambiguous sources. 



in IRAC bands coincides with the ^^CO J = 2 — 1 in- 
tensity. Shocked bright 4.5 /im emission (green) appears 
around most radio continuum sources, except component 
d. 

At least three different velocity components are shown 
in the average spectra of G49.4-0.3 (Figure [5]). The 
brightest component at 52 km s~^ is associated with 
most of the radio continuum sources in G49.4-0.3 be- 
cause th^_jn^aj]__ve]^cit^^f the ionized gas is at 53 km 
s~^ (j Pankonin et al.l [i979f ). However, the velocities of 
the individual components have not been measured. The 
H2CO absorption f or components b and c is observed at 
51 and 63 km s"^ (|Arnal fc Gosslll985 ). and the H I ab- 
sorption for components a, b, c. d, an d f has velocities 
between 42 and 66 km s"^ (jKoolll997[ ). 

Figure [TTl shows the PV map through components a, 

e, and f. In the ^^C0J = 2 — 1 PV map, there are three 
different velocity components. The brightest molecular 
cloud lies at 52 km s~^, another at 62 km s~^ and the 
other at 69 km s~^. All radio sources, except component 

f, are coincident with the molecular clouds at 52 km s~^. 
The H I absorption at 53 km s~^ is dominant toward 
components a, b, c, and d, but for component f it is 
negligible (Kool ll997[ ). Taking into account the dominant 
H I absorption and CO molecular cloud at the velocity 
of 69 km s~^, component f is likely to be associated with 
the molecular cloud at 69 km s~^. 

3.5.3. G49.2-O.3 

Figure [18] shows G49.2-0.3, which is the largest and 
brightest H II region in W51B, in IRAC bands as a 
three-color composite image, ^^CO J = 2 — 1 by thick 
yellow contours, and 21 cm radio continuum by thin 



Fig. 19. — The ^^CO J = 2 - 1 intensity integrated over the 
velocity range from 55 to 75 km s~^ overlaid on the color image 
of the G49.17-0.21 and G49.10-0.27 region composed of IRAC 5.8 
/im (red), 4.5 /im (green), and 3.6 /im (blue). Contour levels are 
20, 40, 60, 80, 100, 120, and 140 K km s~^ Star symbols are 
the radio continuum sources (Koo 1997). YSO candidates (open 
circles) are marked in red for Stage O/I, yellow for Stage II, and 
green for ambiguous sources. 

white contours. The radio continuum emission of G49.2- 
0.3 is apparently coincident with the bright region in 
IRAC bands. The 21 cm continuum brightness decreases 
rapidly toward the molecular cloud. The ^^CO J = 2 — 1 
molecular cloud emission integrated between 57 and 76 
km s~^ encompasses the half of the 21 cm radio contin- 
uum source toward the north. The morphology of the 
radio continuum source and the distribution of the CO 
emission indicate that G49.2-Q .3 is a good example of 
blister type H II region Psraellll978f ). The RRL veloc- 
ity o£_tlie_JO^£egion2s6 7. 6 km s~^ observed in H109Qf 
by iPankonin et al.l (|1979[ ) and the peak velocity of the 
^^C0J = 2 — 1 emission averaged over all pixels shown 
in Figure [T8l is 67.5 km s~^ (Figure [SJ. The H II region 
G49.2-0.3 forms on the boundary of the molecular cloud 
at 67.5 km s~^. The ^^CO J = 2 — 1 emission is strong 
in the east side of the H II region where many YSOs ap- 
pear around {l,b = 49?26, -0?34). The peak velocity of 
^^CO J = 2 — 1 emission is 68 km s~^, which is close to 
the RRL velocity of G49.2-0.3. As discussed in the large 
scale structure of W51B, detailed velocity structure of 
CO lines through G49.2-0.3 and YSO clustering region 
is shown in Figure [9l These two regions divide into two 
different clumps on the ^^CO J = 2 — 1 PV map of Fig- 
ure [9] (mi (i(i/e) and the integrated intensity map of Figure 
[T8l The ratios of ^^/^^i?2-i are smaller than 1.5 at the 
position of H II region and the YSO clustering region 
around 66 km s~^ (Figure [9] bottom). The regions with 
the small ratios are associated with the star-forming ac- 
tivity. 

3.5.4. G49. 17-0.21 and G49. 10-0.27 

G49.17-0.21 and G49.10-0.27 are faint in IRAC bands, 
compared with other radio continuum sources in the 
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Fig. 20. — Spatial and position- velocity diagrams of G49. 17-0.21 
in ^^/^^i?2-i- (a) is channel map at 61 km s"-*^. (b) and (c) 
are position-velocity maps in latitude and longitude across the 
solid lines in (a). Contour levels are 1.5, 3.0, 5.0, and 7.0. (d) 
is composite image of IRAC 5.8 //m (red), 4.5 iim (green), and 
3.6 /im (b lue). The contour s show the 21 cm radio continuum 
emission ( Koo fc MoonI 119971 ). The contour levels are 0.01, 0.03, 
0.05, an d 0.10 Jy b eam~^. The star symbol is the radio continuum 
source (jKod [19971 ). YSO candidates (open circles) on (d) are 
marked in red for Stage O/I, yellow for Stage II, and green for 
ambiguous sources. 

W51B region (Figure [8]). Figure [19] shows the ^^CO 
J = 2 — 1 intensity superposed on the IRAC three-color 
composite image. G49. 17-0.21 is located in the north- 
east of the map. A cometary shape that consists of a 
head around the H II region and a tail toward north- 
east appears in all four IRAC bands as a dark feature. 
The ^^C0J = 2 — 1 integrated intensity map covers the 
bright IRAC emission around the H II regions over a 
large scale. The brightest CO emission coincides with 
the darker part of the IRAC image. In Figure [20l we 
present the ratio of ^^/^^i?2-i on the channel map of 61 
km s~^ and PV maps of G49. 17-0.21. The region in- 
side the first contour (^^/^^i?2-i < 1-5) coincides with 
the dark part of the IRAC image, where YSOs lie. The 
small ratio of -^^/-^^i?2-i implies that the molecular cloud 
around G49. 17-0.21 is exhibiting self-absorption in ^^CO 
J = 2 — 1. In the case of G49. 10-0.27, a dark feature on 
the IRAC image appears around the H II region, where 
many YSOs are identified. The RRL v elocity of 61 km 
s~Ms measured only for G49. 10-0.27 bv iPankonin et aTl 
(IT979I) . 

3.5.5. G49. 1-0.4 

Figure [21] shows the spatial and kinematic structure 
of the molecular cloud around G49. 1-0.4. The molecular 
cloud consists of five clumps on the ^^CO J = 2 — 1 in- 
tegrated intensity map (Figure [2lT a)). Figure [22] shows 
the CO spectra at the peak positions of each clump in 
i^CO J = 2 - 1. Clump 1 has a ^^CO J = 2 - 1 peak 
velocity of ~ 62 km s~^. Clumps 2 and 3 have a ^^CO 
J = 2 — 1 peak velocity of ~ 66 km s~^, which is similar 
to that of the H166a RRL measured from 049. 1-0.4 by 
IPankonin et al.l (|1979[ ). CO emission of clumps 1,2, and 
3 is associated with the dark region on the IRAC compos- 
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Fig. 21. — Spatial and position- velocity diagrams of G49. 1-0.4 
in I^CO J = 2-1. (a) is the ^^CO J = 2-1 intensity integrated 
over the velocity range between 56 and 82 km s"-*^. Contour 
levels are 30, 60, 90, 120, 150, and 180 K km s-^ The grayscale 
is a stretch from to 190 K km s"-*^. Star symbol is the radio 
continuum source (Koo 1997). (b) and (c) are position- velocity 
maps in latitude and longitude across the solid lines in (a). 
Contour levels are 1.6, 3.2, 6.4, 13.0 K. The grayscale is a stretch 
from to 20 K. (d) is the composite image of IRAC 5.8 /im (red), 
4.5 /im (green), and 3.6 /xm (blue). The conto urs show the 21 cm 
radio continuum emission ( Koo fc Moonlfl997I V The contour levels 
are 0.05, 0.10, 0.30, and 0.50 Jy beam~^. YSO candidates (open 
circles) are marked in red for Stage O/I, and green for ambiguous 



ite image and there are several YSOs in the early stage 
of evolution. Near the radio continuum source, there are 
two clumps: clump 4 refers to the 62 km s~^ component 
and clump 5 refers to the 75 km s~^ component. 

It is significant that the -'^^CO J = 2 — 1 inten- 
sity is stronger than the ^^CO J = 1 — in clump 5 
(^^i?2-i/i-o > 2). The high ratio of ^^i?2-i/i-o means 
that clump 5 is significantly different from the other 
clumps. Clump 5 seen in Figure [2TI a) coincides with 
the bright IRAC composite emission and the 21 cm ra- 
dio continuum emission seen in Figure EU^d). The RRL 
velocities in G49. 1-0.4 are measured at 72.4 (H109a), 
70.9 (HlQQg) 71.4 (H137/3) and 68 (H109a) km s"^ 
([Wilson et al.lll970l : IPankonin et"an 119791 : iDownes et all 
ll98Ql ). Assuming the velocity of ionized gas is 72 km s~^, 
clump 5 is associated with G49. 1-0.4 (Figure l2TI b) and 
(c)). Comparing the ratios (^'^^^^R2-i and ^^R2-i/i-o) 
at the ^^CO J = 2 — 1 peak of clumps 3 and 5, the excita- 
tion temperature of clump 5 is about three times higher 
than that of clump 3. The total mass of clump 5 is es- 
timated to be 2.3 X 10^ Mq by a large velocity gradient 
(LVG) analysis using ^^CO and ^^CO J = 2 - 1 lines. 

3.5.6. G49.0-0.3 and G48.9-0.3 

G49.0-0.3 and G48.9-0.3 are located at the west edge of 
the W51B region (Figure [8]). Figure [23] shows the ^^CO 
J = 2 — 1 emission superposed on the IRAC three color 
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Fig. 22.— The ^^CO J = 2-1 (thick sohd hne), ^^CO J = 2-1 
(thin sohd hne), and ■'^^CO J = 1 — (dash-dotted hne) spectra at 
the peak position of each clump in Figure [2]T a). The 75 km s"-*^ 
component in the fourth panel actually belongs to clump 5. 

composite image. IRAC emission is very bright around 
the two radio continuum sources. The bright IRAC emis- 
sion around G49.0-0.3 makes a bow shap e, which has 
been observed as a faint K-band nebulosity (iKumar et al.l 
[2004). Many stellar compone nts seen as white sp ots are 
to the east of the bright bow. iKumar et al.l (|2QQ4[ ) found 
that G48.9-0.3 is composed of two sub-clusters in near- 
infrared J, H, K observations. These two sub-clusters 
also appear in the IRAC three-color composite image. 

The average CO spectra of G49.0-0.3 and G48.9-0.3 in 
Figure [5] have a similar shape. Velocities of G49.0-0.3 and 
G48.9-0.3 at the ^^CO J = 2 - 1 peak intensities are 67.5 
and 68.5 km s~^, respectively. The velocities of the ion- 
ized gas are measured at 64.9 and 64.4 k ni s~^ for G49.Q- 
0.3 a nd G48. 9-0.3 in the H166a line (jPankon in et al." 
[19791 ). The ratios of ^^i?2-i/i-o of G49.0-0.3 and G48.9- 
0.3 are large (> 0.7) around the RRL velocity, which 
indicates that the molecular clouds associated with the 
radio continuum sources are highly excited. 
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Fig. 23. — The ^^CO J = 2-1 intensity integrated over the 
velocity range between 57 and 75 km s~^ overlaid on the color 
image composed of IRAC 5.8 iim (red), 4.5 /im (green), and 3.6 
/xm (blue) around the G49.0-0.3 and 48.9-0.3 region. Contour 
levels are 36, 72, 108, 144, 180, 216, and 25 2 K km s' ^ Star 
symbols mark the radio continuum sources fKod '1997*). YSO 
candidates (open circles) are marked in red for Stage O/I, yellow 
for Stage II, and green for ambiguous sources. 



Many YSOs are found in the relatively dark region of 
the IRAC image along the filamentary structure of the 
molecular cloud. Two bright CO components are coinci- 
dent with two radio continuum sources. We confirm that 
the dense regions are associated with the newly forming 
YSOs by comparing the distribution of YSOs and the 
^^CO J = 2 — 1 intensity on the channel maps. Figure [241 
shows the ^^C0J = 2 — 1 intensity in the velocity range 
from 61 to 75 km s~^ per 1 km s~^. There are three 
groups of YSOs: one to the east of G49.0-0.3, another 
to the west of G48.9-0.3, and the other between G49.0- 
0.3 and G48.9-0.3. Dense molecular clouds associated 
with each group lie in different velocity ranges. YSOs 
between G49.0-0.3 and G48.9-0.3 are coincident with the 
high column density region from the 68 km s~^ channel. 
As seen in the PV diagram of ^^CO J = 2 — 1 in Figured 
(middle), there are three dense molecular clumps around 
G49.0-0.3 and G48. 9-0.3. While the two clumps located 
on either side of the whole cloud are associated with ra- 
dio continuum G49.0-0.3 and G48.9-0.3, the molecular 
clump between the two H II regions lies toward higher 
velocity. These two H II re gions form an e xpanding arc 
of gas in PV diagram (Moon fc Parklll998f ) (Figure I!]). 

4. DISCUSSION 

4.1. Spatial Distribution of the Gas, YSOs, and Radio 
Continuum Sources 

Comparing the velocities of our CO dat a to the veloc- 
ities o f RRL, H2CO and HI absorption (lArnal fc Goss 



'1985; Bieging 1975; "Downes et all |1980r iKa 



Mufson & Liszt 1979; Pan konin et al.l 



19791 ) discussed in 



previous sections, we deduce the spatial distribution of 
the W51 complex. The spatial distribution of the H II 
regions and YSOs around the central active star- forming 
region is shown in Figure [25l Many small clusters of 5-10 
YSOs are apparent. Five of them are marked by large 
circles which happen to be close to PV cuts. Figures [2BI- 
EHlshow the velocity structure of ^^CO J = 2 - 1, ^^CO 
J = 2 - 1, and ^^/^^i^s-i from b = -0?45 to b = -0?20 
by a 0?05 step. We compare the spatial distribution of 
the H II regions and the YSO clusters to the latitude- 
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Fig. 24. Channel maps of G48. 9-0.3 and G49.0-0.3 in the 
■"^^CO J = 2 — 1 emission over the velocity range from 61 to 75 
km s"-*^. Contour levels are 1.3, 2.6, 5.2, and 10.4 K in Tmb- 
Star symbols are the radio continuum sources (Koo 1997). YSO 
candidates (open circles) are marked in red for Stage O/I, yellow 
for Stage II, and green for ambiguous sources. 



velocity (/, v) maps of ^^CO J = 2 - 1, ^^CO J = 2 - 1, 
and their ratio (^^/^^i?2-i) to determine the relative posi- 
tions of the molecular clouds associated with the H II and 
YSO cluster regions. The brightest H II region, G49.5- 
0.4 extends from b = -0?45 to b = -0?30. The mean 
RRL velocity of 59 km s~^ agrees with the strong ^^CO 
and ^^CO J = 2 - 1 intensity. Most of the G49.4-0.3 
sources are located on the b = — 0?30 line; their mean 
RRL velocity is 52 km s~^. As seen on the section of 
b = — 0?30 in Figures [26H28| two strong components 
of molecular gas are located at 52 km s~^ and 68 km 
s~^. Those are associated with G49.4-0.3 and G49.0-0.3. 
The HV molecular components over 70 km s~^ in W51A 
are not evident on the b = — 0?30 line. In Figures [26] - 
[28^ the dense molecular cloud around (l^v) = (49?40, 50 
km s~^) in b = — 0?25 is associated with component of 
G49.4-0.3 and one around {l,v) = (49?20, 61 km s-^) 
inb = -0?20 is connected with G49.17-0.21. The radio 
continuum sources and the associated molecular cloud 
in W51A are behind the HV streaming gas considering 
H2CO absorption line velocities. The HV streaming gas 
in front of G49.5-0.4 is much brighter in CO and more 
extended than that of G49.4-0.3. Spatially, G49.4-0.3 at 
52 km s"^ is located behind G49. 5-0.4 at 58 km s"^ 
Most of the bright molecular clouds in W51B are within 
the velocity range of the HV stream. Toward G49.0-0.3 
and G48.9-0.3, no other molecular clouds are observed 
in the lower velocity range. The sources of W51B within 
the HV streaming gas are located nearest to the Sun. 

YSOs showing clustering are embedded in their par- 
ent molecular cloud. We find the positions of molecu- 
lar clouds associated with YSOs by comparing Figure [25] 
with Figures [26]428] First, many YSOs are concentrated 




Fig. 25. — Spatial distribution of the H II regions (stars) and 
YSOs (small open circles) in the W51 complex, overplotted on a 
color image composed of ■'^^CO J = 2 — 1 intensity integrated over 
velocity range 30-55 km s"-*^ (blue), 56-65 km s"-*^ (green), and 
66-85 km s~-^ (red). The dotted horizontal lines are guide lines 
by 0?05 step. Two vertical dotted lines are guide lines for dividing 
W51A from W51B. Large thick open circles represent the YSO 
clustering regions. 

around {l,b = 49?68, -0?45) in Figure [25] A molecular 
cloud is distributed around {l^v) = (49?68, 70 km s~^) 
in the {l,v) map on b = -0?45 (Figures [MEHl). The 
second YSO cluster around {l,b = 49?62, -0?25) is as- 
sociated with a molecular cloud around (/,v) = (49? 62, 
58 km s-i) on 6 = -0?25 in Figures [26]M] The third 
YSO cluster is located along the b = — 0?35 line between 
G49.5-0.4 and G49.4-0.3 (Figure [25]). In Figure [28] the 
molecular cloud around (l^v) = (49?42, 65 km s~^) at 
b = — 0?35 shows a small ratio of ^^/^^i?2-i where the 
optically thick region coincides with the third YSO clus- 
ter. The fourth YSO cluster is near G49. 2-0.3 at {l,v) 
= (49?28, 65 km s-^) on b = -0?35 in Figures [21^8] 
The fifth YSO cluster is located between G49.0-0.3 and 
G48.9-0.3 in Figure [25] and coincides with the molecular 
cloud at {l,v) = (48?95, 68 km s"^) and b = -0?30 in 
Figures [26]428] A schematic diagram showing the dis- 
tribution of molecular clouds, the H II regions, and the 
YSO clusters is shown in Figure [29] 

4.2. Triggered Star Formation 

In this section, we discuss triggered star formation in 
W51. Triggered star formation was reviewed as three 
distinct mechanisms by Elmegreen (JJ98): (1) direct 
compression of clouds by superno vae, stellar winds, and 
H II regions ([Klein et al.' '1985[), (2) collect and col- 
lapse ([Elmegreen fc La da 1977), and (3) cloud collisions 
(|Loren[[1976 \ 

The arc structure around G49.0-0.3 and G48. 9-0.3 is 
likely to be an example of star formation triggered by the 
collect and collapse process. In a PV diagram through 
the two H II regions (Figure [9]), the molecular cloud com- 
prises a redshifted arc structure with Vlsr = 60 to 73 
km s~^(Moon & Park 1998). Two molecular clumps in 
this arc structure coexist with two H II regions. An- 
other molecular clump is located at higher velocity than 
the locations of the two compact H II regions. Newly 
forming YSOs are associated with the dense molecular 
clumps (Figure [24j) . The shell structure near the W51A 
region (Figure [4j) also shows evidence of star formation 
triggered by the ex pansion of an H II r egion, as discussed 
in an earlier paper (jKang et al.[[2009al ) . Dense molecular 
material has been collected along the shell detected in 
Spitzer IRAC images. YSOs are identified toward the 
densest molecular condensation along the edge of the 
IRAC shell. 
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Fig. 26. — The longitude- velocity map of ^^CO J = 2 - 1 at fixed 
latitudes from b = -0?20 to 6 = -0?45 by 0?05 step. Contour 
levels are 3, 6, 12, and 24 K in Tmb- Two dotted vertical lines are 
guide lines for dividing the W51A and W51B regions. A dotted 
horizontal line shows the boundary of the HV streaming gas. 




Fig. 27. The longitude-velocity map of ^^CO J = 2-1 along 
the latitude from b = -0?20 to 6 = -0?45 by 0?05 step. Contour 
levels are 1.5, 3, 6, and 12 K in Tmb- 
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Fig. 28. — The longitude-velocity map of ^'^/^^R2-i along the 
latitude from b = -0?20 to 6 = -0?45 by 0?05 step. Contour levels 
are 1.5, 3.0, 5.0, and 7.0. 
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Fig. 29. — Schematic diagram of the W51 H II regions, YSO 
clusters, and molecular clouds (compare with the diagram in 
iGo udis & Hippelein (1982)). Dashed circles and filled dots rep- 
resent the H II regions and YSO clusters, respectively. 



iScQville et al.l (| 19861 ) suggested that most high-mass 
OB stars are formed by cloud-cloud collisions based on 
the observational evidence for many GMCs associated 
with H II regions. The studies of star formation trig- 
gered by cloud-cloud collisio ns have been reported in 
otherregions e.g. NGC 1333 (|Lorenlll976l), DR 21-W75 
(jDickel et all goil), W49N (|&r abvn et al. 199 3|), and 
Safflt tarius B2 (jHasegawa et al.i |1994: Mehringer et al.l 
Il993f ). The possibility of triggered star formation 
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by cloud-cloud collisions i n W 5 1 has already been 
suRRested jPa nkonin et afl 119791 : lArnal fc Gossl 119851 : 
iCarpenter fc Sanders 1998). 

Cloud-cloud collisions may compress the interface re- 
gion between two colliding clouds and result in the ini- 
tiation of star formation. The molecular clouds at this 
interface are heated, but the molecular clouds on the 
trailing side remain cold. As a result, self- absorption of 
CO lines is observed toward the interface of cloud-cloud 
collisions. Models for collisions between two different 
clouds have shown that a large cloud can be disrupted 
by the bow shock induced by the collision with a small 
cloud, but the small cloud is compressed by the shock, 
and then th e post-shock g as of the small cloud can form 
a new star (iHabe fc Ohtal ll992f ). 

In Figures [26H281 the velocity of 61 km s~^ divides the 
HV stream from the W51 main components. The HV 
components are related to the streaming motions in the 
Sagittarius spiral arm. Their elongated structure may 
result from the associated spiral-wave shock. We agree 
with earlier suggestions that active massive star forma- 
tion in the W51 H II region complex is due to a collision 
between the molecular clouds in the HV stream a nd the 
molecular cloud below 61 km s~^ (Pankonin et ahl [19791 : 
lArnal fc Go"ssl 119851 : ICarpenter fc~Sanders.,1998E ~In the 
regions G49. 17-0.21 and G49. 5-0.4, the self-absorption 
of ^^CO J = 2 — 1 occurs around 65 km s~^, coincident 
with the dark shadow region in IRAC images. In ad- 
dition, several YSOs are detected inside the low ratio 
regions (^^/^^i?2-i < 1-5). Dark filamentary structures 
in IR bands (IRDCs) are believed to indicate the earli- 
est stages of massive star-forming regions (Beuther et al. 
[2007). The IRDCs associated with G49.17-0.21 and 
G49.5-0.4 have kinematics revealed by the CO maps that 
appear to be consistent with the cloud collision model for 
massive star formation. 

Star formation in the W51 complex appears to have 
occurred in a long and large filament along with the 
H II regions (Figure [1]). We could not find clear evi- 
dence for a time sequence in this formation process from 
our data. YSOs along the filament, however, are younger 
and more massive than those in the surrounding region 
fsee lKang erani2009bL Figures 6 and 7). 

4.3. Star Formation Efficiency and Rate 

Estimating the masses of stars and their parent molec- 
ular clouds is important in understanding how molecular 
clouds evolve and form stars. The star formation effi- 
ciency (SFE) is defined as SFE = M^/{Mgas + M*), 
where Mgas is the mass of cloud a nd M^ is the mass 
of sta rs in the current generation. IWilliams~ fc McKee 
([1997) have found that 10% of the mass of a CMC is 
converted to stars. Koo (1999) reported SFE values of 
7% and 15% in the northwest and southeast components 
ofW51B. 

We can compare the gas mass with the mass currently 
identified as YSOs. We predict that the extended total 
YSO mass of W51 (W51A + W51B) is 2 1 x lO ^Mg, 
using the mass functions from iKang et aTl ()2009b ) by 
extrapolating the fitted power laws down to 1 Mq. Total 
masses of YSOs in W51A and W51B are 0.9 x lO^M© 
and 1.2 X lO^M©, respectively. As discussed in Section 
I we used the LVG formalism for radiative transfer of 



CO lines in molecular clouds to estimate the gas mass. 
Total gas masses of each region are listed in Table [H 
The current YSO mass fractions in W51A and W51B 
are 0.7% and 1.1 %, respectively. 

These estimates are for the mass of YSOs currently 
detected in the selected region by extrapolating in mass 
over our mass function to undetected YSOs. The overall 
star formation efficiency is the ratio of the mass of stars 
formed over the cloud lifetime {Tdoud) to the initial cloud 
mass, or approximately, to the current cloud mass plus 
the mass of all stars formed up to now. We find that only 
about 1% of the cloud mass is currently in YSOs. We as- 
sume the mean timescale for a YSO in the selected mass 
range to be detected by IRAC {Tyso) is ~ 3 x 10^ yr, 
which is the statistical lifetime of high-mas s protostars 
in the Cygnus X molecular cloud complex ( Motte et al.l 
120071 ) . If the final star formation ef ficiency is more like 
10% (| Williams fc McKed [19971 : [K^q 1999) the star for- 
mation process must continue for about 10 x Tyso at 
the present rate to convert 10% of the cloud to stars. So 
the current fraction of YSOs, together with an assumed 
YSO timescale implies the overall time during which the 
cloud is actively forming stars. If the YSO timescale 
for high mass stars is so short, 3 x 10^ yr, for a typical 
cloud lifetime of 3 x 10^ yr and star formation occurred 
at the presently observed rate, the W51 clouds would 
convert 10% of their mass to stars in only 1/lOth of the 
cloud lifetime. We are probably observing these clouds 
at an exceptionally high star formation rate, relative to 
the average over the clou d Hfetime, which typi cally does 
not exceed 5-10 Myrs (jLeisawitz et al.|[l989f ). That is 
consistent with the remark ably high conte nt of massive 
YSOs that we found in Kan g et all (|2009bf ). This "burst 
of star formation" may be a consequence of the passage 
of the W51 CMC through the spiral arm. 

5. SUMMARY 

We observed the molecular clouds associated with the 
active star- forming complex W51 in the ^^CO and ^^CO 
J = 2 — 1 lines with the HHT. We summarize our results, 
combining the new CO data of high sensitivity and res- 
olution, together with Spitzer 3.6 - 24 fim data: 

1. CO emission in the 1?25 x 1?00 area divides into 
three velocity-components, 30 — 55, 56 — 65, and 66 — 85 
km s~^. G49.4-0.3 and some clouds in the northwestern 
part of the map are associated with the lower velocity 
(30 — 55 km s~^) component. The main molecular clouds 
in the velocity range of 56 to 65 km s~^ are distributed 
widely throughout the whole region. The higher velocity 
component (> 66 km s~^) shows an elongated filamen- 
tary structure from southeast to northwest. 

2. The average ratios of ^^i?2-i/i-o and ^^/^^i?2-i 
of all H II regions are greater than the values outside 
the active star- forming regions, which implies that the 
molecular gas directly associated with the H II regions is 
dense and highly excited. 

3. We compare our CO maps with Spitzer dditdi. Strong 
PAH emission near H II regions seen in IRAC bands co- 
incides with bright CO emission, suggesting that those 
molecular clouds are associated with the H II regions. 
Many YSOs are detected in the dark parts near the H II 
regions on IRAC images. 

4. We argue that triggered star formation resulted 
from H II region expansion into clouds in the shell struc- 
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ture near W51A and in the region around G49.0-0.3 
and G48.9-0.3. We also present evidence of star forma- 
tion triggered by cloud-cloud collisions in G49.5-0.4 and 
049. 17-0.21. 

5. We estimate the total gas masses of the various 
interesting regions using an LVO analysis of ^^00 and 
^^00 J = 2 — 1 data. By comparing the total gas mass 
to total mass of YSOs in W51A and W51B, we find cur- 
rent YSO formation efficiencies of 0.7% and 1.1%, re- 
spectively, within a timescale of 3 x 10^ yr during which 
massive stars (> SM©) would be detected as YSOs. The 
current rate of star formation in the W51 OMC should 
then be an order of magnitude higher than the rate av- 
eraged over an assumed cloud lifetime of 3 x 10^ yr, in 
order to convert 10% of the cloud mass to stars. 
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